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The effect of plasma sheaths on the radiation patterns of antennas on hypersonic vehicles is 
theoretically investigated. A simplified mathematical model is employed that consists of an infini- 
tesimal longitudinal magnetic dipole on a plasma-enclosed cylinder. Computed radiation patterns for 
the point source and for a horn and a slotted array arc presented. 

The effects of varying both the permittivity and the loss tangent of the sheath arc considered. 
It is found that in the lossless case the radiated power of an otherwise omnidirectional source is greatly 
reduced in spatial regions beyond the critical angles of the sheath as defined from geometrical optics. 
Also, the transmission loss at broadside due to a lossy sheath is found to be generally greater than 
predicted by plane-wave, plane-sheath approximations. 

These results are extensions of earlier work by Harris. Some of the curves in this paper are 
corrections to the curves in the earlier work. 



1. Introduction 

The presence of ionized sheaths about hypersonic 
re-entry vehicles may have profound effects on the 
characteristics of electromagnetic waves radiated 
from or incident upon such vehicles. Effects such as 
radiation pattern distortion, antenna impedance 
changes, and signal attenuation may result from the 
ionized environment. The magnitudes of these effects 
depend in extremely complicated ways on a combina- 
tion of such physical factors as vehicle shape and 
trajectory, frequency of the electromagnetic waves, 
composition of the atmosphere, etc. It becomes 
important, therefore, to attempt to determine the 
nature of some of these effects by representing the 
actual physical configurations by simplified mathe- 
matical models which are as simple as possible and 
yet can lead to a prediction of the salient features 
of the phenomena of interest. 

This paper is devoted to a study of radiation pattern 
distortion which may result when cylindrical vehicles 
are surrounded by plasma sheaths. Since the primary 
interest is the ability to communicate through the 
plasma, it has been assumed that the radio frequency 
is above the plasma frequency for most cases con- 
sidered. The effects of varying both the relative 
permittivity and loss tangent of the sheath are 
determined. 

The results reported here are an extension of earlier 
work [Harris, 1963]. Radiation patterns in the pres- 
ence of dissipative sheaths and in regions beyond 
the critical angles which result from the presence of a 
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sheath are important additions to this earlier work. 
In addition, the patterns between broadside and the 
critical angle for the lossless plasmas shown in the 
present study are corrections to the patterns included 
in the previous paper. 

2. Presentation of the Problem 

2.1. Mathematical Model 

In this paper an exact solution for the radiation field 
of an antenna in an idealized geometry is presented. 
The geometry is a model for the problem of an antenna 
mounted on the side of a long slender vehicle. De- 
picted in figure 1, the model consists of a finite mag- 
netic line source flush-mounted on an infinitely long, 
perfectly conducting cylinder. The cylinder is sur- 
rounded by a concentric but not contiguous cylindrical 
plasma sheath. This model incorporates some of the 
more salient features of the true problem — the physi- 
cal characteristics of the long slender vehicle, the 
sheath, and the gap between the sheath and the vehicle. 
The low density region beyond the sheath is neglected 
because the propagation frequency is assumed to be 
comparable to the plasma frequency in the high density 
portion of the sheath and therefore much higher than 
the plasma frequency in the low density region. 

The approximation of the actual sheath by a cylin- 
drical sheath with homogeneous properties appears 
to be reasonable, based on approximate solutions for 
electron concentration about re-entry vehicles and on 
the models considered by other investigators [e.g., 
Rotman and Meltz, 1961; Swift, 1964]. The assump- 
tion of infinite lengths, which allows an exact solution 
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FIGURE 1. Geometrical configuration for a line source on a conduct- 
ing cylinder surrounded by a plasma sheath. 



to be obtained, also appears reasonable for the case 
of longitudinal slots, so long as the vehicle and the 
sheath are long compared with their diameters. 

In short, this paper deals with the exact solution of 
a physically meaningful problem. Radiation fields 
for the magnetic line source on the cylinder are com- 
puted as a function of in the plane that contains the 
line source and the cylinder axis. Parameters which 
affect the radiation are the three radii — the cylinder 
radius and the radii of the inner and outer bounderies 
of the sheath (respectively denoted by po, pi, and pi 
in the text) — and the electron density, and collision 
frequency in the sheath. The d-c magnetic effects 
due to the earth's field or any magnetic devices that 
may be in the vehicle are neglected. 

2.2. Selection of Parameter Values 

The relations between the collision and plasma fre- 
quencies and the sheath constitutive parameters are 
[Allis, Buchsbaum, and Bers, 1963] 



and 
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Table 1. Ratio of collision and plasma frequencies to propagation 
frequencies for the sheath constitutive parameters of figures 2 
and 3. 





Fre- 
quency 
ratio 


Relative permittivity. e/e 


0.20 


0.40 


0. 80 


ar/u 6o 
vlto 




0.894 




0.775 


0. 01 
.0167 

.777 


0.05 
.0834 

. 779 


0.10 
.167 
.786 


0.50 

0.834 

1.008 





0. 401 


0.05 
.25 
.414 


0. 10 
.50 
.448 



values of v and a>,, in the sheath for figures 2 and 3 as 
computed from (1) and (2). It should be noted that 
for the most lossy case, cr/co6o = 0.50, the plasma fre- 
quency is greater than the propagation frequency, but 
because of the high collision frequency, the relative 
permittivity, e/6 () = 0.40, remains greater than zero. 2 
The mathematical technique used in determining 
the radiation patterns was taken from Wait [1959] and 
is described in detail in the earlier paper [Harris, 1963]. 
The field expression that yields the radiation field in 
the plane of interest 3 is 



-ikR 



i/co/A) 



R 



X Vi%-h 2 b» 



/>'() COS 



(3) 



where co is the propagation angular frequency, v is 
the collision frequency, and o) }) is the plasma frequency. 
Table 1 has been prepared to indicate the various 



2 Some relationships between flight conditions and parameter values are discussed in the 
report on which this paper is based. [Harris, Villeneuve, and Broca, 1964]. 

3 In the interest of briefness the details are not presented here. For these the reader is 
referred to Harris [1963] and to Harris, Villeneuve. and Broca. [1964]. 
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FIGURE 2. Calculated radiation patterns of an in- 
finitesimal magnetic line source on a cylinder 
enclosed by a lossless plasma sheath with various 
relative permittivities (e/e ) and inner sheath, radii 
(k () pi). 



e.degrees 



a. e/e„ = 0.99. A p, = 3.57r. 

b. e/e„ = 0.90: A p, = 3.5tt. 
-. e/e„ = 0.80; ^p,=47r. 

I. e/e (1 = 0.80; A„p, = 3.5tt. 



e. e/e„ = 0.40; A<,p, = 47r. 

f. e/€o=0.40; A<»p, = 3.5?r. 

g. e/e = 0.20; /t pi = 3.5ir. 
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Figure 3. Calculated radiation patterns of an infinitesimal magnetic line source on a cylinder enclosed by a lossy plasma 
sheath with various relative permittivities (c/eo) and normalized conductivities (cr/co€ ). 

a. e/e = 0.80; o-/cu6o = 0.05. d. e/e = 0.40: o-/we = 0.05. 

b. e/e„ = 0.80; o-/aj€„ = 0.10. e. e/e = 0.40; o-/aje = 0.10. 

c. e/e = 0.40: o-/we„ = 0.01. f. e /e = 0.40; <r/<ue„ = 0.50. 



where R and 6 are the usual spherical coordinates 
measured from the same origin as the cylindrical co- 
ordinates. The effect of the sheath and cylinder is 
entirely represented in the term b m (h) which may be 
written 



b m (h) = 



j: 



dz E. dP 



fm(h). 



(4) 



The integral term in (4) contains the antenna element 
and array factors, while the term fm{h), which is a func- 
tion of the geometry, enters into the cylinder factor. 
This cylinder factor is of principal interest, and is 
discussed in section 3 but is not presented explicitly 
because of its excessive length. The total solution is 
discussed in section 4. 



3. Radiation Patterns 

In this section the computed radiation field in the 
plane of the magnetic line source and the cylinder 
axis is presented. The strength of the source is 
maintained constant for all cases, and the curves 
represent the field of an infinitesimal magnetic dipole 
of fixed dipole moment on the cylinder. Each of the 
curves in figures 2 and 3 may thus be compared 
directly. The field is termed the "cylinder factor" 
and is exclusive of the element and array factors. 

3.1. Lossless Sheath 

A lossless sheath is one in which the conductivity 
is zero. That is, the wave number in the sheath is 
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real and given by 



k — /f(P 



(5) 



indicating that there are no heat losses in the sheath. 
The radiation patterns of a unit magnetic source in 
the presence of such a sheath are shown in figures 
2a through 2g. The dielectric constant in the sheath 
and the sheath thickness are parameters for the curves. 
The critical angle C shown on the patterns is defined 
by 



sin 2 r = €/€<> 



(6) 



The values kopo = 27r and kop2 = 6ir are constant and 
correspond, respectively, to a cylinder radius of 1 free 
space wavelength and an outer sheath radius of 3 free 
space wavelengths. The values of 3.57T and \ir for 
/c pi, the inner sheath radius, correspond respectively 
to 1.75 and 2 free space wavelengths. 

The e/e { )= 0.99 case (fig. 2a) is essentially the sheath- 
free case and exhibits none of the amplitude fluctua- 
tions which are characteristic of the radiation in the 
presence of denser sheaths. The case e/eo— 0.90 
(fig. 2b) corresponds to a ratio of plasma frequency to 
propagation frequency of 0.316 and is somewhat more 
characteristic of sheath effects. Three fairly modest 
peaks appear between 10 and 22°, or in the vicinity 
of the critical angle which is at 18° in this case. The 
power drops to more than 20 dB below broadside power 
at about 15°, excluding the peak at 11°, while in the 
e/€o=0.99 case it drops to this cutoff value at 9°. 

The €/€ = 0.80 cases are shown in figures 2c and 
2d. These correspond to sheath thicknesses {p2~ P\) 
of 1 and 1.25 free space wavelengths, respectively. 
The latter case has rather sharp cutoff at about 22.5° 
and a rather interesting narrow beam at 13.5° of approx- 
imate width 1° and magnitude 26 dB. Such a beam is 
accounted for physically as the radiation field of a 
leaky wave which, due to leakage of energy as it propa- 
gates, is attenuated slowly along the cylindrical struc- 
ture. 4 The pattern for the one-wavelength-thiek sheath 
exhibits deeper nulls than that for the 1.25-wavelength 
sheath. While the pattern essentially remains above 
the —20 dB level from 90 to about 20.5°, signal levels 
fall below the — 20 dB level over 1° widths at angles as 
great as 33°. 

The e/eo — 0.4 cases (figs. 2e and 2f) exhibit a great 
number of peaks and valleys. Of interest here, in 
addition to the narrow beams, is the great difference in 
the patterns caused by the 1/4-wavelength change in 
sheath thickness. Even the level at broadside varies 
by 4 dB in the two cases, and while the deviations from 
a smooth curve occur at about the same angles, these 
deviations generally differ in form in the two cases. 
The pattern essentially drops below —20 dB at 40.5° 
in the 1.25-wavelength case and at 42.5° in the 1 -wave- 
length case. 



The final lossless case presented corresponds to 
e/€o = 0.2 (fig. 2g). Here the pattern essentially drops 
below —20 dB at 61°, and once again a spectrum of 
narrow, high-intensity beams is present. It is evident 
that decreasing the sheath permittivity increases the 
number of significant leaky wave beams. 



3.2. Lossy Sheath 



a. Pattern Characteristics 

A lossy sheath is one in which the conductivity is 
finite and nonzero. In this case the propagation con- 
stant in the sheath is given by 



: hi 



€ 



. (T 
O)€o 



(7) 



4 The relationship hctween leaky wave, i.e., nonspectral, modes associated with a struc- 
ture and the characteristics of the radiation patterns from the structure is discussed in 
great detail in Tamir and Oliner [1962]. 



which implies that there is electromagnetic energy 
lost to the sheath in the form of heat due to collisions 
between the electrons and other particles in the sheath. 
The assumption made here is that the power absorp- 
tion does not raise the temperature of the sheath suffi- 
ciently to alter the electron density. Figures 3a 
through 3f indicate the radiation patterns of an infini- 
tesimal magnetic dipole on the cylinder in the presence 
of a sheath of thickness 1.25 free space wave- 
lengths with relative permittivities e/e<) = 0.80 and 
e/eo = 0.40 and with normalized conductivity cr/coeo 
as a parameter. 

Figures 3a through 3f, taken in conjunction with 
figures 2d and 2f, indicate a phenomenon that had 
been previously hypothesized: that the sharp peaks 
of the pattern in the lossless case would be flattened 
out as a result of increased conductivity. This pre- 
dicted smoothing is apparent in the patterns. Indeed, 
the radiation field in the highest loss cases is virtually 
monotonic. 

b. Power Loss 

The radiation patterns shown in this section cannot 
be used directly to determine the power radiated by 
the antenna nor the power loss in the sheath since 
such calculations would require a knowledge of the 
antenna impedance [Villeneuve, 1965]. If the assump- 
tion is made, however, that the impedance seen by 
the antenna is essentially constant for fixed sheath 
thickness and permittivity, then the difference in 
the areas under the power curves for lossless and lossy 
sheaths represents the power lost in the sheath. Of 
perhaps greater interest is the reduction in transmitted 
power at broadside as the losses are increased, again 
under the assumption that the computed patterns for 
a given sheath thickness and permittivity represent 
the same constant power out of the antenna. Table 
2 indicates this power loss at broadside for the cases 
of figures 3a through 3f. Also included in this table 
for comparison is the power density transmitted 
through a plane sheet when a plane wave is normally 
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Table 2. Broadside transmission loss through plasma sheath from the following expression [Stratton, 1941]: 



Relative permittivity, e/e« 


Normalized conductivity. cr/a>€o 





0.01 


0.05 


0.10 


0.50 


0.80 
A* 
B* 

0.40 
A* 
B* 











-0.80 
- .30 


-2.0 
-2.4 

-3.6 
-0.67 


-4.0 
-3.9 

-6.4 
-3.7 


-24.8 
-22.6 



*A = ratio of power transmitted in dB at broadside through a lossy cylindrical sheath to the 
power transmitted through a lossless cylindrical sheath of the same thickness (1.25 
free space wavelengths) and the same permittivity. 

*B = ratio of power density transmitted through a plane lossy sheath to the power trans- 
mitted through a plane lossless sheath of the same thickness (1.25 free space wave- 
lengths) and permittivity when a plane wave is normally incident on the plane sheath. 

incident. The thickness and the constitutive param- 
eters of the plane sheath correspond to those same 
values in the cylinder problem. Powers are shown 
relative to the lossless sheath. It is intersting to note 
that the approximation of a plane wave propagating 
through a plane sheath, which is often used to estimate 
transmission loss, gives fairly comparable results 
with the transmission loss at broadside for the cylinder 
case and thus is somewhat justified as a first order 
approximation, for the cases considered here although 
the plane approximation is seen to generally yield a 
smaller loss than does the cylindrical geometry. 
Different conclusions have been reached by Swift 
[1964] for the overdense plasma and smaller values 
of kopo. The plane wave approximation was computed 



^ e -ik (f(k/k -l) 



& in 



2 + f + ^ 
/Co k . 



k-ko \ 2 
k+k 0i 



(8) 



where k is defined in (7), k is the free space propaga- 
tion constant, and d is the thickness of the sheath. 



3.3. Comparison With Planar Radiation 
Problem 

As one method of checking on the reasonableness 
of the radiation patterns for the cylindrical case, the 
limiting case of a slot in a groundplane was considered. 
A plasma slab was located above the groundplane 
and separated from it as illustrated in figure 4. On 
applying the boundary conditions at the groundplane 
and at the slab surfaces and on using the principle of 
stationary phase [Di Francia, 1955] to evaluate the radi- 
ation field, the following expression is obtained for 
E^ in the plane (/> = 0. 



= 4 cos 2 6 



Ve 







(9) 



where 6 is measured from broadside and e 2 is the slab 
relative permittivity. The denominator is given by 



I e ifco(*2-*i>vV s in2 [(co s 2 6 - sin 2 6 + e 2 + 2cos 6 Ve 2 - sin 2 fl)e**o *i cos 9 + (1 - € 2 )e-* fc o*i cos e ] 
\- c-*tf*2-*i> V€ 2 -sIn2 9 ^ cog2 0_ sin 2 + €2 -2cos dVZ - sin 2 6)e ik o?i cose + (l ~ e 2 )e- ik o*i cos ] 



(10) 



TO FIELD POINT 




GROUNDPLANE 



When sin 2 is greater than e 2 , the square root is nega- 
tive imaginary. This expression for E. has been com- 
puted as a function of 6 for the case of e 2 = 0.20, 
&oZi = 1.57T, and koZo = 477. This case corresponds to 
the cylindrical case for which &oPo = 2?t, kopi = 3. 577, 
and kop2 = 6tt. The resulting curve normalized to its 
value at = is shown in figure 5. While this curve 
is much smoother than those for the corresponding 
cylindrical geometry of figure 2g, several very im- 
portant characteristics are evident. First, most of 
the radiation is confined between broadside and the 
critical angle. Second, the peak which occurs near 
the critical angle corresponds to that type which occurs 
when there is no gap between the slab and the ground- 
plane. [Tamir and Oliner, 1962]. Third, a high, very 
narrow peak occurs at about = 57.9°. This peak 
results from the presence of the interface at z\ and dis- 
appears as z\ approaches zero. The value of the field 
at this angle, denoted by O , is given by 



\ E 4 (0o)| =cos 6o e fe o<*2-*i) Vsin 2 o -e 2 . 



(ID 



FIGURE 4. Slot in plasma-covered groundplane. 
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Figure 5. Radiation pattern of slot in plasma-covered groundplane. 



It may be shown that, as Z\ approaches zero, do ap- 
proaches the critical angle, given by 



sin 2 6 C = e 2 . 
Then the expression for E^Oo) reduces to 



|£*(0o)| = V(l-€,) 



(12) 



(13) 



and the peak has disappeared. 

The feature of importance in the radiation pattern 
is the existence of the very narrow peak which results 
from the interface at z\. It tends to confirm the exist- 
ence of these peaks for the cylindrical geometry where 
the analytical expressions are much more difficult to 
interpret. The existence of many more peaks in the 
cylindrical case is attributed to the more complicated 
nature of waves propagating in that configuration. 



4. Discussion of Results 
4.1. Leaky Wave Effects 

The radiation patterns of figures 2 and 3 have two 
principal characteristics: increasing numbers of nar- 
row peaks and valleys with decreasing permittivity, 
and suppression of the peaks with increasing conduc- 
tivity. Both characteristics may be explained in 
terms of the presence of leaky waves [Tamir and 
Oliner, 1962]. 

Leaky waves can be described physically as electro- 
magnetic waves which propagate along a structure 
while losing power to space at a uniform rate. Each 
of the m modes in (3) has its own set of leaky wave 
poles and, since one point source generates all of the 
modes with equal amplitude, the curves of figures 2 
and 3 represent a direct superposition of the fields 
contributed by the leaky waves as well as what might 
be described as the direct radiation from each of the 
modes. The fact that each of the peaks represents a 



leaky wave field of a single mode is evidenced by term- 
by-term computer print-outs of the summation in (3). 
In the cylinder problem at hand, a determination 
of the leaky wave propagation constants involves the 
solution of an extremely complicated transcendental 
equation. This solution has not been attempted, but 
the computed field patterns indicate, on one hand, 
that a lowering of the sheath permittivity increases 
the number of significant leaky wave poles and, on the 
other hand, that the pole positions are sensitive to 
sheath thickness since quarter-wave variations in 
sheath thickness for the e/eo = 0.8 and e/eo — 0.4 cases 
(figs. 2c, d, e, and f) yield significant pattern changes. 
The leaky wave concept affords no simplification of the 
radiation problem for the lossless case because a de- 
termination of the leaky wave pole positions is a more 
difficult task than actual computation of the radiation 
patterns. When the sheath is lossy, however, the 
leaky wave concept may be useful since the general 
form of the suppression of peaks may be derived by 
considering the leaky wave to deliver power to the 
sheath as well as to the far field. Thus, the first order 
effect of nonzero conductivity in the sheath is to alter 
the z component of the wave number from the form 
[^-h/'-nJ to [j>+/(t7 + 8)] where 8, which is the measure 
of the power loss to the sheath, is primarily dependent 
on the conductivity. This increase in the imaginary 
part of the propagation constant tends to yield a 
smaller, more flattened beam radiated by the leaky wave 
[Harris, 1963]. When trie loss is sufficiently great, 
beams no longer appear and the radiation field is 
monotonic. This result is evident in the computed 
patterns for the cases of largest sheath conductivity. 

4.2. Pattern Considerations 

As has been previously stated, the radiation pat- 
tern of an antenna on a space vehicle is the product 
of the cylinder factor, element factor, and array factor. 
The curves of figures 2 and 3 are cylinder factors 
for various sheath configurations and constitutive 
parameters. If it is desired to determine the radiation 
patterns of antennas in the indicated environments, 
then it is sufficient to add the antenna pattern, speci- 
fied in dB, to the various cylinder factors shown in 
this report. 

If the general results of such a pattern addition are 
considered, the most important factors are: For high 
altitude flight, when the sheath losses are small, 
there will be a blackout region from approximately 
the critical angle to endfire. A blackout region is 
defined as one in which the cylinder factor is generally 
down by 20 dB from broadside although spurious 
beams with broadside power will be present in the 
blackout region. A narrow-beam broadside array 
will thus operate satisfactorily from the pattern point 
of view, 5 while an omnidirectional antenna will be 
cutoff at the critical angle. An endfire array would 
not be expected to be constructed of slots parallel to 
the cylinder axis. 



"'This statement excludes impedance considerations; see Villeneuve [1965] for example. 



1341 



£ 20 



30 



















(a) 










































e c =CR 

AN 

I 


ITICAL 
GLE 
1 













\ 


—N 


I 


to 










(b) 








V 


1 


1 ' 
















1 






i 
















U 


i 


V 


1 










6 C CRITICAL 
















ANGLE 

1 i 








A. 





90 8 ° 7 ° 60 50 40 30 20 10 90 80 70 60 50 40 30 20 10 

©.degrees 

FIGURE 6. Radiation patterns of wide-beam horn. 

a. Measured pattern of horn alone. 
b. Calculated pattern of horn on cylinder enclosed by a lossless plasma sheath. e/e ( ) = 0.4, <x = 0, kopo = 2TT, koPi = S.5ir, knp-> = 6Tr. 
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FIGURE 7. Radiation patterns of 20 -element, 25-dB Tchebyscheff array. 

a. Measured pattern of array alone, 
b. Calculated pattern of array on a cylinder enclosed by a lossless plasma sheath. e/e = 0.4, cr = 0, kop„ = 2Tr. fcopi = 3.57T, k p-> = 67r. 
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As an example, patterns for a horn antenna and a 
20-element array are shown in figures 6 and 7; actual 
patterns measured in the absence of the cylinder and 
plasma are presented in figures 6a and 7a, respec- 
tively, while the radiation pattern determined analyti- 
cally for each type of antenna by multiplying the meas- 
ured patterns by the cylinder factors is shown in 
figures 6b and 7b. The horn and the narrow-beam 
array are assumed to be on a vehicle for which e/e () 
= 0.4, cr = 0, and A'opi — 3.57T as in figure 2f. The 
effect of the cylinder factor in the total patterns can 
be seen. The pattern of the narrow-beam array re- 
mains relatively unaffected by the plasma while that 
of the wide-beam horn is significantly tapered. (Only 
points at quarter-degree intervals were used in the 
computation of these composite patterns.) 

When a vehicle is at intermediate altitudes, the loss 
factor due to increased collision frequency in the 
sheath plays an important role. Various loss factors 
for typical sheath conditions are shown in the rows 
marked "A" in table 2. As has been pointed out, 
the plane-wave plane-sheath approximation gives a 
general idea of the transmission loss through the sheath 
at broadside. 



5. Conclusions 



The principal results of this 
antenna patterns presented in 
They indicate the effect of 
radiation patterns of antennas 
by an illustration of these effec 
matical model. An important 
establishment of a rule-of-th 
factor falls off rapidly beyo 



study are the numerous 

the body of the report, 
plasma sheaths on the 

on hypersonic vehicles 
ts on a simplified mathe- 

additional result is the 
umb that the cylinder 
nd the critical angles. 



The transmission loss through a lossy sheath at broad- 
side is generally greater than the result predicted 
by the usual approximation of a plane wave at normal 
incidence on a plane sheath, although this approxi- 
mation is useful in obtaining an estimate of the loss. 
The computations for this report (table 2) show a 
maximum difference in transmission loss of 3 dB be- 
tween the cylindrical and the plane-wave values. 
Finally, there is the large number of narrow beams 
due to leaky wave poles which are generated by the 
sheath. 
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